J. Pure Appl. & Ind. Phys. Vol.2 (3), 222-233 (2012) 



A Semi-Circle Theorem in Couple-Stress Fluid in the 
Presence of Magnetic Field in a Porous Medium 

AJAIB S. BANYAL 

Department of Mathematics, 
Govt. College Nadaun, Hamirpur-177033, H. P., INDIA. 

(Received on : May 17, 2012) 
ABSTRACT 

A layer of couple-stress fluid heated from below in a porous 
medium is considered in the presence of uniform vertical 
magnetic field. Following the linearized stability theory and 
normal mode analysis, the paper through mathematical analysis of 
the governing equations of couple-stress fluid convection with a 
uniform vertical magnetic field in porous medium, for any 
combination of perfectly conducting free and rigid boundaries of 
infinite horizontal extension at the top and bottom of the fluid, 
established that the complex growth rate <7 of oscillatory 
perturbations, neutral or unstable for all wave numbers, must lie 

inside a semi -circle |cr| = £ 2 Q 2 , in the right half of a complex 

<7 -plane, where Q is the Chandrasekhar number and£ is the 
porosity of the porous medium, which prescribes the upper limits 
to the complex growth rate of arbitrary oscillatory motions of 
growing amplitude in the couple-stress fluid heated from below in 
the presence of uniform vertical magnetic field in a porous 
medium. The result is important since the exact solutions of the 
problem investigated in closed form, are not obtainable for any 
arbitrary combinations of perfectly conducting dynamically free 
and rigid boundaries. The result is also in accordance with 
corresponding configuration of couple-stress fluid when the 
porosity £ = 1 , Banyal 1 . 

Keywords: Thermal convection, Couple-Stress Fluid, Magnetic 
field, PES, Chandrasekhar number. 
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1. INTRODUCTION 

Stability of a dynamical system is 
closest to real life, in the sense that 
realization of a dynamical system depends 
upon its stability. Right from the 
conceptualizations of turbulence, instability 
of fluid flows is being regarded at its root. A 
detailed account of the theoretical and 
experimental study of the onset of thermal 
instability (Benard Convection) in 
Newtonian fluids, under varying 
assumptions of hydrodynamics and 
hydromagnetics, has been given by 
Chandrasekhar 2 and the Boussinesq 
approximation has been used throughout, 
which states that the density changes are 
disregarded in all other terms in the equation 
of motion, except in the external force term. 
The formation and derivation of the basic 
equations of a layer of fluid heated from 
below in a porous medium, using the 
Boussinesq approximation, has been given 
in a treatise by Joseph 3 . When a fluid 
permeates through an isotropic and 
homogeneous porous medium, the gross 
effect is represented by Darcy's law. The 
study of layer of fluid heated from below in 
porous media is motivated both theoretically 
and by its practical applications in 
engineering. Among the applications in 
engineering disciplines one can name the 
food processing industry, the chemical 
processing industry, solidification, and the 
centrifugal casting of metals. The 
development of geothermal power resources 
has increased general interest in the 
properties of convection in a porous 
medium. Stommel and Fedorov 4 and 
Linden 5 have remarked that the length scales 
characteristic of double -diffusive convecting 



layers in the ocean may be sufficiently large 
so that the Earth's rotation might be 
important in their formation. Moreover, the 
rotation of the Earth distorts the boundaries 
of a hexagonal convection cell in a fluid 
through porous medium, and this distortion 
plays an important role in the extraction of 
energy in geothermal regions. The forced 
convection in a fluid saturated porous 
medium channel has been studied by Nield 
et al. 6 . An extensive and updated account of 
convection in porous media has been given 
by Nield and Bejan 7 . 

The effect of a magnetic field on the 
stability of such a flow is of interest in 
geophysics, particularly in the study of the 
earth's core, where the earth's mantle, which 
consist of conducting fluid, behaves like a 
porous medium that can become 
conductively unstable as result of 
differential diffusion. Another application of 
the results of flow through a porous medium 
in the presence of magnetic field is in the 
study of the stability of convective 
geothermal flow. A good account of the 
effect of rotation and magnetic field on the 
layer of fluid heated from below has been 
given in a treatise by Chandrasekhar 2 . 

MHD finds vital applications in 
MHD generators, MHD flow-meters and 
pumps for pumping liquid metals in 
metallurgy, geophysics, MHD couplers and 
bearings, and physiological processes such 
magnetic therapy. With the growing 
importance of non-Newtonian fluids in 
modern technology and industries, 
investigations of such fluids are desirable. 
The presence of small amounts of additives 
in a lubricant can improve bearing 
performance by increasing the lubricant 
viscosity and thus producing an increase in 
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the load capacity. These additives in a 
lubricant also reduce the coefficient of 
friction and increase the temperature range 
in which the bearing can operate. 

Darcy's law governs the flow of a 
Newtonian fluid through an isotropic and 
homogeneous porous medium. However, to 
be mathematically compatible and 
physically consistent with the Navier-Stokes 
equations, Brinkman 8 heuristically proposed 

the introduction of the term —V 2 q, (now 

£ 

known as Brinkman term) in addition to the 



Darcian term — 



V^l J 



q . But the main effect 



is through the Darcian term; Brinkman term 
contributes very little effect for flow through 
a porous medium. Therefore, Darcy's law is 
proposed heuristically to govern the flow of 
this non-Newtonian couple-stress fluid 
through porous medium. A number of 
theories of the micro continuum have been 
postulated and applied (Stokes 9 ; Lai et al. 10 ; 
Walicka 11 ). The theory due to Stokes 9 allows 
for polar effects such as the presence of 
couple stresses and body couples. Stokes's 9 
theory has been applied to the study of some 
simple lubrication problems (see e.g. Sinha 



et al} 2 ; Bujurke and Jayaraman 13 ; 



Lin 14 ). 



According to the theory of Stokes , couple- 
stresses are found to appear in noticeable 
magnitudes in fluids with very large 
molecules. Since the long chain hyaluronic 
acid molecules are found as additives in 
synovial fluid, Walicki and Walicka 15 
modeled synovial fluid as couple stress fluid 
in human joints. The study is motivated by a 
model of synovial fluid. The synovial fluid 
is natural lubricant of joints of the 
vertebrates. The detailed description of the 



joints lubrication has very important 
practical implications; practically all 
diseases of joints are caused by or connected 
with a malfunction of the lubrication. The 
external efficiency of the physiological joint 
lubrication is caused by more mechanisms. 
The synovial fluid is caused by the content 
of the hyaluronic acid, a fluid of high 
viscosity, near to a gel. A layer of such fluid 
heated from below in a porous medium 
under the action of magnetic field and 
rotation may find applications in 
physiological processes. MHD finds 
applications in physiological processes such 
as magnetic therapy; rotation and heating 
may find applications in physiotherapy. The 
use of magnetic field is being made for the 
clinical purposes in detection and cure of 
certain diseases with the help of magnetic 
field devices. 

Sharma and Thakur 16 have studied 
the thermal convection in couple-stress fluid 
in porous medium in hydromagnetics. 
Sharma and Sharma 17 have studied the 
couple-stress fluid heated from below in 
porous medium. Kumar and Kumar 18 have 
studied the combined effect of dust particles, 
magnetic field and rotation on couple-stress 
fluid heated from below and for the case of 
stationary convection, found that dust 
particles have destabilizing effect on the 
system, where as the rotation is found to 
have stabilizing effect on the system, 
however couple-stress and magnetic field 
are found to have both stabilizing and 
destabilizing effects under certain 
conditions. Sunil et al. 19 have studied the 
global stability for thermal convection in a 
couple-stress fluid heated from below and 
found couple-stress fluids are thermally 
more stable than the ordinary viscous fluids. 
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Pellow and Southwell 20 proved the 
validity of PES for the classical Rayleigh- 
Benard convection problem. Banerjee 
et al. 21 gave a new scheme for combining the 
governing equations of thermohaline 
convection, which is shown to lead to the 
bounds for the complex growth rate of the 
arbitrary oscillatory perturbations, neutral or 
unstable for all combinations of dynamically 
rigid or free boundaries and, Banerjee and 
Banerjee 22 established a criterion on 
characterization of non-oscillatory motions 
in hydrodynamics which was further 
extended by Gupta et al. 23 . However no such 
result existed for non-Newtonian fluid 
configurations, in general and for couple- 
stress fluid configurations, in particular. 
Banyal 24 have characterized the non- 
oscillatory motions in couple-stress fluid. 

Keeping in mind the importance of 
couple-stress fluids and magnetic field in 
porous media, as stated above,, the present 
paper is an attempt to prescribe the upper 
limits to the complex growth rate of 
arbitrary oscillatory motions of growing 
amplitude, in a layer of incompressible 
couple-stress fluid in a porous medium 
heated from below, in the presence of 
uniform vertical magnetic field, opposite to 
force field of gravity, when the bounding 
surfaces are of infinite horizontal extension, 
at the top and bottom of the fluid and are 
perfectly conducting with any combination 
of dynamically free and rigid boundaries. 
The result is important since the exact 
solutions of the problem investigated in 
closed form, are not obtainable, for any 
arbitrary combination of perfectly 
conducting dynamically free and rigid 
boundaries. 



2. FORMULATION OF THE PROBLEM 
AND PERTURBATION EQUATIONS 

Here we consider an infinite, 
horizontal, incompressible electrically 
conducting couple-stress fluid layer, of 
thickness d, heated from below so that, the 
temperature and density at the bottom 
surface z = 0 are T n and 0 n and at the 



upper surface z = d are T d and p d 
respectively, and that a uniform adverse 



temperature gradient f3 



dT 



dz 



is 



maintained. The fluid is acted upon by a 
uniform vertical magnetic field 

H{0,0,H). This fluid layer is flowing 
through an isotropic and homogeneous 
porous medium of porosity £ and of medium 
permeability k A . 

Let p, p, T, 77 , // c and q{u,v, w) 

denote respectively the fluid density, 
pressure, temperature, resistivity, magnetic 
permeability and filter velocity of the fluid, 
respectively Then the momentum balance, 
mass balance, and energy balance equation 
of couple-stress fluid and Maxwell's 
equations through porous medium, 
governing the flow of couple-stress fluid in 
the presence of uniform vertical magnetic 
field (Stokes 9 ; Joseph 3 ; Chandrasekhar 2 ) are 
given by 
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The equation of state is 

p = p 0 [l-a{T-T 0 )], (6) 

Where the suffix zero refer to the values at 





= -V 




+ 8 








KPoj 




V PoJ 



Po 



\7tp 0 



(l) 



(2) 



AT -> 

E — + ( q .V)T = *y 2 T, (3) 



V.H=0. 



(4) 



£^ = (H.V)q+£T ? V 2 H, (5) 



d__(^ -i v 

Where , — + & Q • » , stands for 

dt dt 

the convective derivatives. 



Here 



E = £ + {\-£) 



\Po C v J 



is a constant, while P s , C s and P 0 , C v , 

stands for the density and heat capacity of 
the solid (porous matrix) material and the 
fluid, respectively, £ is the medium 

— > 

porosity and T (x, y, z) . 



the reference level z = 0. Here g(0,0,—g) is 
acceleration due to gravity and (X is the 
coefficient of thermal expansion. In writing 
the equation (1), we made use of the 
Boussinesq approximation, which states that 
the density variations are ignored in all 
terms in the equation of motion except the 
external force term. The kinematic viscosity 

V , couple-stress viscosity jU , magnetic 

permeability jU e , thermal diffusivity K , and 

electrical resistivity TJ , and the coefficient of 

thermal expansion CC are all assumed to be 
constants. 

The basic motionless solution is 

q = (0,0,0) , p = p 0 Q. + oflz), 

p=p(z), T = —f3z + T Q , (7) 

Here we use the linearized stability 
theory and the normal mode analysis 
method. Assume small perturbations around 
the basic solution, and let dp , dp , 6 ' , 

q(u,v,w) and h = (h x ,h y ,h z ) denote 

respectively the perturbations in density p , 

— > 

pressure p, temperature T, velocity q (0,0,0) 

and the magnetic field H = (0,0, H ) . The 
change in density dp , caused mainly by the 
perturbation 0 in temperature, is given by 
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p + Sp = p o [l-a{T + 0-T o )] 
= p-ap Q 0 

i.e. dp = -ap o 0 . (8) 



Then the linearized perturbation 
equations of the couple-sress fluid reduces to 



1 d_q_ 
dt 

Me 



—VSp-lae— 

£ dt p a Kj 



Vx/z \xH 



(9) 



Where k x ,k are the wave numbers along 
the x- and y-directions, respectively, 

k = {k x 2 + ky 2 ) 2 , is the resultant wave 

number, n is the growth rate which is, in 
general, a complex constant and, 
W(z),&(z) and^(z) are the functions of z 
only. 

Using (14), equations (9)-(13), 
Within the framework of Boussinesq 
approximations, in the non-dimensional 
form transform to 



V.<? = 0, 

E — = ]8w + i& 2 6> 
dt 



v./Uo. 



d_h 
dt 



f 



H.V 



(10) 



(ii) 



(12) 



q+£J]V 2 h_ ( i3) 



(D 2 -a 2 



£ + P, 



(D 2 -a 2 ) 



i j 1 1 



W = -Ra 2 e + QD(D 2 -a 2 )K 

(d 2 -a 2 - p 2 ct)k = -DW , 
and 

(d 2 -a 2 -Ep 1 <j)® = -W , 



(15) 



(16) 



(17) 



Where we have introduced new 
coordinates [x',y',z') = (x/d, y/d, z/d) in 
new units of length d and D = d I dz' ■ For 
convenience, the dashes are dropped 
hereafter. Also we have substituted 



a = kd,<7 



lid V . 

, /?,= — , is the thermal 

V K 



3. NORMAL MODE ANALYSIS 

Analyzing the disturbances into 
two-dimensional waves, and considering 

disturbances characterized by a particular Prandtl number; p 2 =—, is the magnetic 
wave number, we assume that the 
Perturbation quantities are of the form 



[w,0,h z \=[w{z),©{z),K(z)]exp 

(ik x x + ik y y + nt), (14) 



Prandtl number; 



> d 2 



is the 



dimensionless medium permeability, 

jU/(p 0 d 2 ) 
b = , is the dimensionless 
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couple-stress viscosity 



parameter; 



gaBd 4 

R = , is the thermal Rayleigh 



KV 



2 j2 



number and Q 



M e H z d 



is the 



Chandrasekhar number. Also we have 
Substituted 



( /?j2 \ 



W=W m , 0 



J3d 



v K J 



<h£> 



Km, 



a.ndD @ =dD, and dropped (©) for 
convenience. 

Now consider the case for any 
combination of the horizontal boundaries as, 
rigid-rigid or rigid-free or free-rigid or free- 
free at z=0 and z=l, as the case may be, and 
are perfectly conducting. The boundaries are 
maintained at constant temperature, thus the 
perturbations in the temperature are zero at 
the boundaries. The appropriate boundary 
conditions with respect to which equations 
(15)-(17), must possess a solution are 

W = 0= © , on both the horizontal 

boundaries, (18) 

DW=0, on a rigid boundary, (19) 
D z W=0, on a dynamically free 



boundary, 



(20) 



K = 0, on both the boundaries as the 
regions outside the fluid are 
perfectly conducting, (21) 

Equations (15)-(17) and 

appropriately adequate boundary conditions 



from (18)-(21), pose an eigenvalue problem 
for £7 and we wish to Characterize <T ; , 

when <7 r > 0 . 

4. MATHEMATICAL ANALYSIS 

We prove the following theorems: 

Theorem 1: If R> 0, F>0, 00, (7 r > 0 
and <T i ^ 0 then the necessary condition for 
the existence of non-trivial solution 
{W,®,K) of equations (15) - (17) and the 
boundary conditions (18), (21) and any 
combination of (19) and (20) is that 

\<j\(eQ . 

Proof: Multiplying equation (15) by W* 
(the complex conjugate of W) throughout 
and integrating 

the resulting equation over the vertical 
range of z, we get 



— + — 1\w'{d 2 -a 2 )wdz-— ]w'{d 2 -a 2 fwdz 
y £ ^i Jo Pi o 

1 1 

-Ra 2 jWQdz + Q \w*d{d 2 -a 2 )Kdz > ( 22 ) 



F V. 



Taking complex conjugate on both sides of 
equation (17), we get 

(d 2 - a 2 - E Pl a )®* =-W\ (23) 

Therefore, using (23), we get 
i i 

jWQdz = - \®{d 2 -a 2 -E Pl a )&*dz > ( 24 ) 

0 0 

Also taking complex conjugate on both sides 
of equation (16), we get 
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[d 2 -a 2 —p <j*]k* =-DW* (25) Therefore, using (25) and using boundary 

2 condition (18), we get 



}w*d(d 2 -a 2 )Kdz = -\DW*(d 2 -a 2 )Kdz = } K~(d 2 - a 2 \D 2 - a 2 - p 2 a*)K*dz , (26) 

0 0 0 



Substituting (24) and (26) in the right hand side of equation (22), we get 

— + — IjV^D 2 -a 2 )tfdz-— \\V*(d 2 -a 2 ) 2 Wdz 

l l 
= Ra 2 \Q(D 2 - a 2 - Ep l a*^) t dz + Q \k(d 2 - a 2 \D 2 - a 2 - p 2 a*)K*dz , (27) 



Integrating the terms on both sides of equation (29) for an appropriate number of times by 
making use of the appropriate boundary conditions (18) - (21), we get 



\{DW\ 2 +a 2 \w\ 2 \lz + — \(\D 2 W\ 2 +2a 2 \DW\ 2 +a 4 |wf \fe 



= Ra 2 j{/)(-) (-) -/./> a (-) '}/- 

o 

D K + 2a 2 DK -a : K )l:-(Jp.(J \(dK ' + a 2 K ')lz 



(28) 



And equating the imaginary parts on both side of Eq. (28), and cancelling <J i 0) 
throughout, we get 

1 \{dw\ 2 +a 2 \w\ 2 \lz = -Ra 2 Ep l j\&\ 2 dz + Qp 2 \(l)K + a 2 K ')/:■ (29) 



We first note that since W and K satisfy W(0) = 0 = W(l) and ^(0) = 0 = ^(1) 
in addition to satisfying to governing equations and hence we have from the Rayleigh-Ritz 
inequality 25 
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i i 

j\DW\ 2 dz>x 2 \\w\ 2 dz, (30) 



and 



\\DK\ 2 dz>7T 2 \\K\ 2 dz, 



(31) 



Further, multiplying equation (16) and its complex conjugate (25), and integrating by parts 
each term on right hand side of the resulting equation for an appropriate number of times and 
making use of boundary conditions on K namely K(0) = 0 = ^(1) along with appropriate 
boundary condition (18) - (21), we get 

i i , , i i 

J](D 2 -a 2 )K\ 2 dz + 2p 2 a r \\DK\ 2 + a 2 \K\ 2 )dz+ p 2 \(J\ 2 \\K\ 2 dz = \\DW\ 2 dz , (32) 



since <J r > 0 , O l ^0 therefore the equation (32) gives, 



J](D 2 -a 2 )K\ 2 dz(j\DW\ 2 dz, 



(33) 



o 

And 

i 



(34) 



\\K\ 2 dz{—^\\DW\ 2 dz, 

It is easily seen upon using the boundary conditions (21) that 

)\DK\ 2 +a 2 \K\ 2 )dz = Real part of j_ j^Z) 2 -a 2 )Kdz\ ^ )k*(d 2 - a 2 )Kdz 

o { o Jo 

< \k*(d 2 -a 2 )K\lz < \k*\\(d 2 -a 2 )K\iz, 

0 0 

1 1 
= j\K\\(D 2 -a 2 )K\dz <|j|z| 2 Jz| 2 ||(D 2 -a 2 )z| 2 dz| 2 , (35) 

(Utilizing Cauchy-Schwartz-inequality) 
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Upon utilizing the inequality (33) and (34), inequality (35) gives 
\{dk\ 2 +a 2 \K\ 2 ]dz<-^- r \\DW\ 2 dz' 

Now R ) 0 and Q)0 , utilizing the inequality (36), the equation (29) gives, 



(36) 



Q_ 
\o\ 



\\DW\ 2 dz + ^- \\w\ 2 dz + Ra 2 E Pl \\®\ 2 dz <0, 



2 1 



(37) 



and therefore , we must have 



\a\(eQ. 
Hence, if 

O r > 0 and £7, * 0 , then \o\{eQ . 

And this completes the proof of the theorem. 

5. CONCLUSIONS 

The inequality (40) for O r > 0 and 
O i ■£ 0, can be written as a 2 + a 2 {e 2 Q 2 , 

The essential content of the theorem, from 
the point of view of linear stability theory is 
that for the configuration of couple-stress 
fluid of infinite horizontal extension heated 
form below, having top and bottom 
bounding surfaces are of infinite horizontal 
extension, at the top and bottom of the fluid 
and are perfectly conducting with any 
arbitrary combination of dynamically free 
and rigid boundaries, in the presence of 
uniform vertical magnetic field parallel to 
the force field of gravity, the complex 
growth rate of an arbitrary oscillatory 
motions of growing amplitude, lies inside a 
semi -circle in the right half of the <T r <T i - 



(38) 



plane whose centre is at the origin and radius 
is equal to£Q , where Q is the 
Chandrasekhar number and£ is the porosity 
of the porous medium. The result is 
important since the exact solutions of the 
problem investigated in closed form, are not 
obtainable, for any arbitrary combinations of 
perfectly conducting dynamically free and 
rigid boundaries. 
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